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Abstract 



o 

We study primordial density perturbations generated by the late decay of a curvaton field whose 
decay rate may be modulated by the local value of another isocurvature field, analogous to models 



of modulated reheating at the end of inflation. We calculate the primordial density perturbation 
and its local-type non-Gaussianity using the sudden-decay approximation for the curvaton field, 
recovering standard curvaton and modulated reheating results as limiting cases. We verify the 
Suyama-Yamaguchi inequality between bispectrum and trispectrum parameters for the primordial 
density field generated by multiple field fluctuations, and find conditions for the bound to be 
saturated. 
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I. INTRODUCTION 



Inflation in the very early universe provides a classical cosmology that drives the uni- 
verse towards spatial flatness and homogeneity It also provides an origin for primordial 
density perturbations through the quantum fluctuations of light scalar fields, stretched by 
the inflationary expansion to super- Hubble scales. Originally structure was assumed to orig- 
inate from fluctuations in the inflaton field driving inflation, but more recently it has been 
realised that there are many possible mechanisms through which scalar field fluctuations 
could generate the observed primordial density perturbations Examples incude the late 
decay (some time after inflation has ended) of a curvaton field modulated reheating 

at the end of inflation, where the rate of inflaton decay is modulated by the local vacuum 
expectation value (VEV) of an isocurvature field ]6, 71, or an inhomogeneous end of infla- 

nn 

tion [8|, |9[. All of these mechanisms for the origin of the primordial density perturbations 
may be distinguished from adiabatic perturbations in the inflaton field driving inflation as 
they introduce local- type primordial non-Gaussianity [10], most simply characterised by the 



non-linearity parameter, f NL 
becomes small 12]. 



Ill ] . which typically becomes large when the transfer efficiency 



In this paper we go beyond the simplest curvaton models to include the possible modu- 
lation of the curvaton decay by a modulator field, Xi m analogy with modulated reheating 
at the end of inflation. In modulated reheating scenarios the rate of decay of a massive field 
a is assumed to be dependent on the VEV of a modulator field, x- Consider a Lagrangian 
which describes the decay of the cx-field oscillations into ^-particles 

C = - V {a) - \(x)<tW> - U( X ) . (1) 

If the field, a, is displaced from its minimum during inflation, then the field oscillates once 
H < m a where m 2 = V"(a) at its minimum. Note that the potential, V(a), may deviate 
from a simple quadratic potential, but we will assume it can be described by a simple 
quadratic potential once it begins to oscillate about its minimum, with effective mass m a . If 
the effective mass of the x field remains small while the a field oscillates, |i7"(x)| < H 2 < m 2 , 
then the decay rate of massive cr-particles into light ^-particles is given by T oc A 2 which is 
assumed to be a function of the local VEV of x- F° r example, if the coupling A is a linear 
function of x then the decay rate is a quadratic function T oc x 2 - 
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Any light field during inflation acquires an almost scale-invariant distribution of pertur- 
bations due to small-scale vacuum fluctuations being stretched up and frozen-in on super- 
Hubble scales. In particular, we will assume that light fields acquire a Gaussian distribution 
of perturbations at Hubble-exit (k = (aH)*) with power spectrum 



Local variations in the x field change the local rate of reheating and hence the primordial 
density perturbations, which can be represented by (, the metric perturbation on uniform- 
density hypersurfaces [lsj in the primordial radiation-dominated era, some time after infla- 
tion has ended. 

In the original modulated reheating scenario a is the inflaton whose oscillations dominate 
the energy density of the universe immediately after inflation comes to an end. In this work 
we will consider the case where a is a curvaton field whose energy density is sub-dominant 
during inflation. We assume its effective mass is small compared to the Hubble scale dur- 
ing inflation and hence it also acquires an almost scale-invariant spectrum of fluctuations. 
The curvaton mass becomes larger than the Hubble scale either immediately at the end 
of inflation or some time after, when H = m a and the curvaton field eventually decays. 
Fluctuations in the local VEV of the curvaton during inflation lead to perturbations in the 
amplitude of oscillations and the density of curvaton particles when the field oscillates after 
inflation. These density perturbations are transferred to the radiation when the curvaton 
decays. 

If both the modulator field and the curvaton field are light during inflation then their 
vacuum fluctuations can generate a primordial density perturbation, in addition to any 
adiabatic density perturbations produced from adiabatic fluctuations in the inflaton field, 
</>, during inflation. 

The outline of this paper is as follows. In section [III we describe how we can estimate the 
primordial density perturbations produced by the modulated curvaton decay in terms of the 
inhomogeneous densities and metric on an instantaneous decay hypersurface, generalising 
previous results for non-linear perturbations from curvaton decay with an unmodulated de- 
cay rate [14]. We calculate observable quantities including the primordial power spectrum, 
tensor-scalar ratio, bispectrum and trispectrum. We show that the Suyama-Yamaguchi in- 
equality between the tree-level bispectrum and trispectrum holds and consider the condition 




(2) 
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for this inequality to be saturated. We present our conclusions and discussion in section IHIl 



II. SUDDEN-DECAY APPROXIMATION 

We will work in the sudden-decay approximation where the curvaton decay is modelled 
as an instantaneous transfer of energy from the curvaton field oscillations, p CT , into radiation, 
p.y. This has been shown to be a good approximation to the full numerical results in the 

n 

usual curvaton scenario |14| . 

If the curvaton decay rate is a constant, T = f , then the decay hypersurface, H dcc = T, 
corresponds to a uniform-density hypersurface (on super-Hubble scales at the decay time) 
with 

p dec = 3M P1 2 f 2 . (3) 

In the presence of fluctuations in the local VEV of x, the modulator field, and therefore 
local fluctuations of the decay rate, we have 

p dcc = 3M P1 2 r 2 ( X ) • (4) 

We will allow for the existence of inhomogeneous perturbations of the radiation density, 
C 7 , after inflation but before the curvaton decays, due to adiabatic inflaton field fluctuations 
during inflation, C7 = Cci = —H5<f>/<j). Before the curvaton decay, the curvature perturbation 
on uniform-radiation-density hyper surf aces, £ 7 , and uniform-curvaton-density hypersurfaces, 



( a , are independently conserved [l5|. After the decay, assuming the decay products are rela- 



tivists, the curvature perturbation on uniform-total-density hypersurfaces, (, is conserved. 
On the decay hypersurface itself, we therefore have Q, Il6 ] 

Pcr,dcc Pcr,dccC ^ ^ , P7,dec P^,dec& ^ 7 ^ ^ j Pdcc Pdec^ ^ j (5) 

where (dec is the curvature perturbation (5N from a flat hypersurface) on the decay hyper- 
surface. Note that we neglect the energy density of the modulator field \ throughout. 
Matching p 7 + p a = p at the decay time yields 

(1 - fi CT , dec )e 4 ^-^ + ^dece 3 ^"^ = e 4 «"^ , (6) 



where 



n a = ^. (7) 
p 
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A. Linear perturbations 



At linear order, one can expand the above relation ([6]) to give 

- (dec ^ TZ~ [3ft CT ,dccCx + 4 ( X - ^,dec)C 7 ~ 4 C] (8) 
" '<r,dec 

Using the above relation to eliminate (dec in © an d also using 01]), results in 

C*Cy-/^ + /(Cr-Cr), (9) 
where we have defined the transfer parameter 

/^T%^- (10) 

The change in the curvature perturbation, C — C7 i n Eq. due to the modulated curvaton 
decay is seen to arise from the relative entropy perturbation between the curvaton density 
and the radiation density, ( a — £ 7 , and the perturbed decay rate, 5T ~ T'Sx- In the case of 
homogeneous curvaton decay rate, 5T = 0, we reproduce the standard curvaton result. The 
inhomogeneous decay rate adds an extra term in the primordial curvature perturbation (J9]) 
proportional to /. 

Note that since x remains overdamped throughout the decay, we assume that we can 
neglect its background evolution and hence its perturbation on the decay hypersurface cor- 
respond directly to its perturbation on spatially flat hypersurfaces during inflation. The 
curvaton, a, is overdamped during inflation, but starts to oscillate when H = m a . We 
assume the curvaton density is still negligible at this time, so the curvaton field fluctuations 
on this surface correspond to curvaton density perturbations on a uniform radiation-density 
hypersurface 

p a = j9 ff e 3 «*-« . (11) 

Allowing for possible evolution of the curvaton field from the end of inflation up until the 
point at which it starts oscillating we write 

Pu = l -mlg 2 {a + 5a). (12) 



At linear order we thus have 



Sa = 3(( a -Q = 2 9 -^. (13) 



where a + 5a describes the local curvaton VEV at the end of inflation. 
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B. Non-linear perturbations 



Expanding (jlj) and (jSJ) to second order yields 

r 



c-c 



dec 



ilnf£ 



2 v r / ~ 2F Sx + 4 I r r 2 ' 5x 



r" 



Eliminating £d ec m ® and solving for ( order by order yields 



./ 



r 



C ~ C 7 + - f —6 X 



/(i-/)(3 + /) o2 r/(i-/)(3 + /) 



lc 



+ 



./ 



6r 

3^ + ^(9-/(/ + 2)))5 X 2 



-5 



isr 



(14) 



(15) 



36 V I" ' 2T 2 
where we define the entropy isocurvature perturbation 



ST 



S a = 3(C CT - C 7 ) - 2-^ - - 



9' 



1 



l9_ 

n l2 



bo 1 



(16) 



Observables 



When calculating observables, such as the power spectrum and higher-order correlators 
of the primordial density perturbations, it will be convenient to express the non-linear per- 
turbation ( in terms of the perturbed logarithmic expansion, N = j Hdt, from an initial 
spatially flat hypersurface where the scalar field perturbations originate during inflation and 



3,Q 



a final uniform- density hypersurface during the subsequent radiation-dominated era 16 
we note that the total curvature perturbation can be rewritten by 

C = NM a + \N ab 5<p a H b + \N abc 5<p a H b H c + -. (17) 
where 5<ft a are the three Gaussian fields 5<ft, 5\ and 5a , N a = N a = dN/d<ft a and N ab = 

N ab = d 2 N/d(j) a (i) b . 

Comparing Eq. (fT7|) with ( 1T5]) we identify 



2M P1 % 



3<? 



6r 



A' 



/(i-/)( 3 + /)rv 



iV™ = -4- 





9r<7 
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r g"9 
1 + — - 


4/ 


2/ 2 




3 
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3 




/ 
36 


r'T 

9 6 


-2/- 


-/ 2 





(18) 

(19) 
(20) 
(21) 



in which 



^-(fn " ( 22 ) 



Note that higher-derivatives with respect to the inflaton field <fi can be neglected during 
slow-roll inflation. 

From Eqs. (ll8HT9]) we identify 

f,4 = ^ , (23) 
da 3 g 

_ 9/ /(i-/)( 3 + /) r 
/x -%" 6 r' 1 j 



This will allow us to calculate all higher derivatives of N starting from Eqs. f fl8lfT9|) . One 
can verify this is consistent with Eqs. ( 1201I2TT) . 

The power spectrum is given, at leading order, by 

n = n, + t v „ +/ > (£) V ^ (I + 1 + £) (|); (25) 

in which, in analogy with the inflaton contribution, we have defined [3] 

-4 (75k) 2 • ^ i8 (/fk) 2 (26) 



The relative contribution of each field to the power spectrum (1251) is given by the weights 
w a defined via V^ a = w a V( in which 



N 2 e- 1 
Wa= N} + W + N* = ^ + tf + e?- (27) 



Note that w$ + w a + w x = 1. 
The spectral index is then 

n s - 1 = w x (n Cx - 1) + w ff (n CCT - 1) + (1 - w x - w a )(n u - 1) (28) 

The tensor to scalar ratio is also given by 

r = 16e x w x = 16e a w a = 16(1 — w a — io x )e^ . (29) 

The local- type primordial bispectrum is characterised at leading order in term of /jyx 

6 N N b N ab 

5 Jnl = jkw (30) 
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FIG. 1: A logarithmic plot for /tvl as a function of relative energy density /. In both plots we 
assumed g oc a and neglected the contribution from inflaton field {w^ ~ 0). For the left plot we 
also assumed T oc x 2 anci in the right one, we have set w a = 0.2. The apparent singularity in the 
left plot is due to a change of sign in /nl- 



which we can write as 

fNL = wlfffL + 2w X W ^NL + W lfNL , 

where we write the different contributions as 



6 ra 

~^JNL 



fO~X. — 
JNL — 
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N x N a 
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x 



1 
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i 
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3fl + £?)-2/-/- 



2 V 9 
1-/K3 + /) 
/ 

2 r'T 

9,1 -3I^ 



(31) 

(32) 
(33) 
(34) 



where we have used Eqs. (118H2ip . 

In the limit % — > 1 then /jvl -> f^ L and we recover the standard result for the curvaton 



19]. In the opposite limit where w x — > 1 and / — > 1 we recover the standard result 



for modulated reheating 



1 - 



r'T 

p2 



The primordial trispectrum is composed of distinct two terms 
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9nl 
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54 (N a N a ) 



(35) 

(36) 
(37) 
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FIG. 2: A plot for TjVL as a function of relative energy density /. The values of the independent 
parameters are the same as in FigfT] 



which in our case we can write as 
25 



36 



+w a w x (w a + w x ){f N \f + 2w a wlt/ L f* L + wKf^L? 
9nl = wlg a NL + 3wlw x g a N a * + 3w a w 2 x g°™ + w x g* L 

where we identify the different contributions to g^i as 
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4/ 2 V 9 



9 fg"'g 2 



,9" 9 
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/ 



g"g 



— (^7r + 3 ^]--( 1 + ^]+^( 1 - 9 ^)+ 1 °/ + 3/ 
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J9 
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25 
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25 5jVi 



AT. 



AT3 > 



r'T 



p///p2 



135 - 54/ - 22 f + 10/ d + 3/ 4 - 18(9 - 2f - f 2 )— + 36 



(38) 
(39) 



(40) 



(41) 



(42) 



(43) 



t2 \ J v J J > p/2 1 p/3 

where we have used Eqs. (l20H24p . 

In the curvaton limit w a — > 1 then c/tvl — ^ fi^L and we recover the standard result for the 
curvaton jljj]. In the opposite limit, w x — >■ 1, we recover the result for modulated reheating 
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FIG. 3: A plot for g^i as a function of relative energy density /. Again the values of the indepen- 
dent parameters are the same as in FigQ] 



when / 1 (lOj, |22| 



x 50 
9nl ~ g 



2-3— + 



p2 



F /3 



(44) 



D. The Suyama-Yamaguchi inequality 



Here we verify that the Suyama-Yamaguchi (SY) inequality 22J, stating that t^l > 
|/atl) 2 , holds at tree-level in our model x . 
A direct analysis shows that 

25 



in which 



Cl 



<(l-w x ) 



c 2 = 2WH(l-2™ x )/£* 



(45) 

(46) 
(47) 



c 3 = u£(l - + w a w x (w a + w x - ±w„w x )f3* + 2w x w 2 a (l - 2 W(T )r NL f N x L (48) 

We wish to determine the sign of K in Eq. ([43]) which is a quadratic function of f^ L for 
ci 7^ 0. Hence we re- write Eq. (145|) as 



« c 2 ^ A 



where 



* = «.|/Si+jf) +-, 



A = whu a {\ -w a - w x )(w a fx L + w x f a N x L f 



(49) 



(50) 



1 Note that loop-corrections violate this tree-level inequality 23, [24 1 
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Both the coefficient, c\ in Eq. ( )46|) . and the discriminator, A in Eq. ( )50|) are non-negative, 
since the weights, w a and w x , and their sum, w a + w x , are bounded between zero and one. 
Hence K given in Eq. is non-negative and we conclude that t nl > (|/atl) 2 as required. 

One may ask under what conditions the SY inequality is saturated. Firstly this can 
occur if ci = and C2f NL + C3 = 0, which requires either w a = 1 or w x = 1, i.e., either 
the curvaton fluctuations or the modulator fluctuations dominate the primordial density 
perturbations corresponding to effectively a single source for the primordial density field, or 
w$ = 1 corresponding to the Gaussian case, t^l = Inl = 

For ci 7^ 0, the SY equality is saturated when A = and 2ci/^ L + c 2 = in Eq. 
This either requires 

wlfNL = KfNL = -^xInI (51) 

which implies Tnl = Inl = 0, or requires w a + w x = 1 and 

Wo Tnl ~ W xflVL = ( W <r ~ ™ X )fNL ( 52 ) 

or, equivalently, 

fX f°X fa _ f°X 

_ TNL TNL _ JNL TNL (rn\ 

fa 1 XX _0/°"X ' /a , /X _9/°"X ' 

•/JVL iTVL NL JNL JNL ^J NL 

,sib 

greater than f^ L , i.e. 



Note that this is possible only when f^ L and f^ L are either both less than or both 



iflVL ~ fSL) ifjSTL ~ fNl) > • (54) 



III. DISCUSSION 



The curvaton and modulated reheating scenarios have previously been studied as distinct 
models for the origin of structure from quantum field fluctuations during inflation. Here we 
have considered a curvaton scenario where the curvaton decay rate may be modulated by a 
second scalar field. Thus fluctuations in two independent fields may be responsible for both 
the primordial power spectrum and primordial non-Gaussianity described by higher-order 
correlations. 

The relative contributions to the primordial power spectrum ( 1251) are determined by 
the weights, w a and w x defined in (1271) . The overall contribution of both curvaton and 
modulator fluctuations to the first-order primordial density perturbation is proportional to 

11 



the fractional density in the curvaton at the time of decay, / defined in ({TO]) . The relative 
contribution depends on the fractional perturbations in the curvaton amplitude of oscillation, 
g, versus the decay rate, T: 

^-fyf- (55) 

We recover (i) previous results for the curvaton scenario in the limit that curvaton field 
fluctuations dominate the primordial power spectrum (w a ^> w x ) and (ii) previous results 
for modulated reheating in the limit that the curvaton dominates the total energy when it 
decays (/ — > 1) and modulated fluctuations dominate (w x ^> w a ) . 

We also allow for an adiabatic density perturbation produced due to inflaton perturba- 
tions during inflation, Q,. This is required to be Gaussian, hence the inflaton contributes to 
the first-order density perturbation, but not to higher-order correlators. Any detection of 
local-type primordial non-Gaussianity would be evidence of non-adiabatic field fluctuations 



25). 



playing a role in the origin of large-scale structure, see however 

Non-Gaussianity can become very large, either for / <C 1 or due to non-linear evolution, 
g"g ^> g' 2 , in the curvaton scenario. Indeed it is bounded by current observations j^J 

f NL = 37.2 ± 19.9 . (56) 

Such a large value is difficult to achieve in standard modulated reheating unless the decay 
rate is strongly dependent upon the modulator field, T"T /T' 2 ~ 8. On the other hand one 
can easily have large non-Gaussianity in a modulated curvaton reheating when / <C 1 even 
when the modulator dominates the primordial fluctuations, u> x — >- 1. 

In the simplest scenario of a quadratic curvaton potential, leading to linear evolution of 
the curvaton field, g" = 0, and a linear coupling, A(x) in Eq. (JTJ, leading to a quadratic 
decay rate, T oc % 2 , the primordial non-Gaussianity is a function of /, w a and w x . Hence a 
measurement of the three lowest-order non-Gaussian correlators, /jvx, tnl and gNL, would 
be required to determine /, w a and w x . 

There are some general predictions for the simplest model of a quadratic curvaton po- 
tential with g" = 0. The contributions to the primordial bispectrum, given in Eq. (|3T|) . 
from fff£ and /jyx, defined in Eqs. (!3~4]l and ( |33l) . are non- negative for any T"T < T' 2 , which 
includes a quadratic decay rate. Thus in the modulated curvaton decay when g" = 0, we 
find a lower bound /nl > /jvx — — 5/4, generalising the result found previously for a single 
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curvaton [ijj], and multiple curvaton decays 27]. This bound is saturated, Jnl = — 5/4, for 
w a = 1 and / = 1. 

For this simple quadratic curvaton scenario with linear evolution, such that g" and g'" can 
be neglected, we find the third-order trispectrum parameter g a NL oc f^ 1 oc /jvl whereas the 
second-order trispectrum parameter t nl oc f~ 2 oc f^ L . Hence t^l is much larger than g a NL 
if /jvl is large. On the other hand g™* oc f~ 2 oc f% L even for g" ~ 0, and hence gNb/^NL is 
not necessarily suppressed for a simple quadratic curvaton in a modulated curvaton scenario 
with 6T + 0. 

We generally expect g NL oc r^x oc /| rf due to non-linear evolution of the curvaton field 



with a self- interaction potential 2814321] . On the other hand the self-interacting curvaton can 



give rise to strongly scale-dependent non-Gaussianity [331 ]. while the modulated curvaton 
decay with a quadratic curvaton potential gives rise to non-Gaussianity which is scale- 
independent. 

We have verified that the trispectrum parameter t^l obeys the Suyama-Yamaguchi in- 

n, 

equality [22| tnl > (36/25) f% L . It is saturated when the curvaton perturbations (w a —> 1) 
or modulator perturbations (w x — > 1) dominate, or in the trivial case of Gaussian perturba- 
tions when — > 1. It can also be satisfied for particular parameter values even when both 
curvaton and modulator fluctuations contribute to the primordial density field, w a + w x = 1 
given by Eq. f )53|) . This emphasizes that a single-source for the primordial density field is 
a sufficient condition for the SY inequality to be saturated, but not a necessary condition 
22]. 

Throughout this work we have assumed that the decay products of the curvaton rapidly 
thermalise leaving no residual isocurvature perturbations. If the matter asymmetry inherits 
a different density perturbation from the overall radiation density (because it comes ex- 
clusively from either the curvaton decay products or the pre-existing radiation before the 
curvaton decay, but not both) then it may leave a residual matter isocurvature perturba- 



tion 



34 



ll-l36l|. This would be further evidence of the origin of structure from non-adiabatic 
field perturbations during inflation, and measurements of the relative amplitude of residual 
isocurvature perturbations and their correlations with the adiabatic density perturbation 
could give independent constraints on model parameters. It would be interesting to investi- 
gate whether these could then distinguish modulated curvaton decay from standard curvaton 
or modulated reheating scenarios. 
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FIG. 4: Plots for the Suyama-Yamaguchi equality /inequality (j45]). In both figures Wa, = 0. Left: 
/ = 0.2 with generic value of fft L , f$ L and ffi L . Right: The particular case ffr L = f$ L = 5f^ x L is 
plotted for arbitrary w a . As expected the SY inequality is saturated when w a = 0, 1 or when Eq. 
is met with w a + w x = 1, which in this particular case, yields w a = w x = 0.5. 



Although we have assumed that the modulator field has negligible energy density, it 
would be interesting to consider possible observational signatures of the eventual decay of 
ihe modulator field, assuming that it too eventually decays into standard model particles 



37j |. analogous to our recent study of the effect of the late-time decay of a field responsible 



for the inhomogeneous end of hybrid inflation |38j] . We leave this for future work. 

Note added: While completing this work, we became aware of related work by Langlois 



and Takahashi 



Both papers should appear on the arXiv on the same day. 
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